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The iron-chromium aUoy and its derivatives are widely used for their remarkable resistance to cor- 
rosion, which only occurs in a narrow concentration range around 9 to 13 atomic percent chromium. 
Although known to be due to chromium enrichment of a few atoms thick layer at the surfaces, 
the understanding of its complex atomistic origin has been a remaining challenge. We report an 
investigation of the thermodynamics of such surfaces at the atomic scale by means of Monte Carlo 
simulations. We use a Hamiltonian which provides a parameterization of previous ab initio results 
and successfully describes the alloy's unusual thermodynamics. We report a strong enrichment in 
Cr of the surfaces for low bulk concentrations, with a narrow optimum around 12 atomic percent 
chromium, beyond which the surface composition decreases drastically. This behavior is explained 
by a synergy between (i) the complex phase separation in the bulk alloy, (ii) local phase transitions 
that tune the layers closest to the surface to an iron-rich state and inhibit the bulk phase separation 
in this region, and (iii) its compensation by a strong and non-linear enrichment in Cr of the next few 
layers. Implications with respect to the design of prospective nanomaterials are briefly discussed. 

PACS numbers: 64.75.Nx, 68.35.bd, 61.66.Dk, 64.70.kd, 81.30.Bx 



The iron- chromium alloy and its derivatives are in- 
expensive, have satisfactory mechanical properties and 
above all exhibit a remarkable resistance to corrosion: it 
is the most widely used class of alloy in the world. Its out- 
standing corrosion resistance is known for a century^ to 
only occur in a narrow range of concentrations, around 10 
atomic percent of chromium (at. % Crj2!. Their excellent 
properties make the m c andidate materials for future fu- 
sion nuclear reactor^^lil^ one of the reasons that induced a 
considerable amount of work on the various aspects of the 
Fe-Cr alloy both experimentall}- ^ ^ and theoretically^''. 

Corrosion resistance of stainless steels is due to the pas- 
sivation of the material by an inert, chromium rich layer 
at the interface between the alloy and the environment, 
i.e. at the surfaces. Passivation is a phenomena inher- 
ent to how much Cr is located at the surfaces, which 
is a non-linear function of the bulk concentration-. In 
austenitic Fe-Cr, which only exists at high temperatures 
above ~ 800 C and for less than w 10 at. % Cr, the more 
chromium in the bulk, the more chromium in the surface 
and thus the more stainless the alloy. In ferritic Fe-Cr 
alloys, the picture is more complex. Without additive 
elements, the Cr content at which the alloy is passivated 
is narrow, from 9 to 13 at. % Cr, beyond which occurs 
an increase in the corrosion rate and a strong decrease in 
mechanical properties. 

This important property of stainless steels has been 
extensively studied, but its complex origin at the atomic 
scale has remained a missing understanding, subject to 
controversial findings: How chromium causes passiva- 
tion, i.e. how it interacts and reacts with chemical el- 
ements coming from the environment like dioxygen or 
hydrogeiP^ is out of the scope of this study. The reader 
is refered to Greeley et for a review of surface chem- 
istry of metal surfaces at the atomic and electronic scale. 
Surface reaction requires nevertheless that Cr is present 



in large enough quantity on the surface to form a few 
atoms thick protective layer, e.g. of chromium(III) oxide 
Cr203. How chromium atoms enrich the surfaces remains 
unclear. Venus and HeinriclP^ shew by angle-resolved 
Auger electron spectroscopy that Cr atoms deposited on 
a whisker of Fe (100) migrates from the surface to the 
first few layers, in contradiction with the expected ten- 
dency. This surface-alloying has been clearly identified to 
be linked to anomalies in the magnetic properties of the 
Cr/Fe system, specifically the change in surface magneti- 
zation at low Cr cover age an d the strong interactions be- 
tween surface Cr atomJI^MI, Ropo et al. showed by First 
Principles that a pure Fe-Cr surface behaves like stain- 
less steels with respect to Cr enrichmenlj^. They also put 
in evidence a competition between the relative stabilities 
of the surfaces and the complex thermodynamics of the 
bulk alloy. Later, ah initio calculations revealed that un- 
expected interactions between subsurface Cr atoms and 
surface Fe atomJ^ are at the origin of an anomalouJ^^'^ 
segregation behavior of Cr in Fe in the dilute regime. 

At temperatures of industrial and technological inter- 
est, i.e. between 300 and 600 K, the body-centered cu- 
bic (bcc) solid-solution of Fe-Cr shows a miscibility gap 
from 9 - 13 to 94 - 99 at. % Cr^. Inside, a phase- 
separation occurs into an iron-rich bcc solid solution, a, 
and chromium-rich bcc precipitates, q', as one would ex- 
pect for a binary alloy that seemed to have a segrega- 
tion tendency, i.e. that mix solely for entropic reasons. 
However, both theoretical and experimental studies sub- 
verted this simple picture, showing favorable dissolution 
energy of chromium in iron up to an anomalously high 
« 7 at. % CjEsHSSIdT^e to a competition between repulsive 
Cr-Cr interactions and attractive Fe-Cr interactiond^iltMl, 
The increased chromium content leads to more frustrated 
magnetic interactions in between Cr atoms that make 
the dissolution exothermic at low concentrations, then 
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endothermic. Several theoretical model^^'^lHlll have suc- 
cessfully reproduced the sign change of the mixing energy. 

Here, the iron-chromium ferritic stainless steel is mod- 
eled by the Hamiltonian proposed recently in Ref.'^Sl. It is 
specifically designed to reproduce both (i) the whole ex- 
perimental a-a' phase diagram at all temperatures and 
compositions, and (ii) the change of sign of the mixing en- 
ergies. It is also compatible with large-scaled simulations, 
because of its conceptual simplicity and its rigid bcc lat- 
tice nature: The internal energy AHnux — —^Cb (1 — Cf,) 
is a function of the bulk concentration in chromium, ct,, 
and of the order energy ^l described in terms of local- 
concentration and temperature dependent pair interac- 
tions: 

^-E^i^'^ + ^s-^e^), (1) 

i 

where z^*-* is the coordination number of shell i and ej^l is 
the pair interaction between atoms of type j and j' on ith 
neighbor sites. Influenced by the strategy of Caro et al^, 
the order energy is advantageously expressed as a sim- 
ple concentration and temperature dependent Redlich- 
KisteJ2fi! expansion: 

nix,T)^{x^a){l3x^+jx + S)(^l-^y (2) 

where x is the local concentration and T the temperature. 
Discrete mixing energies have been calculated ab initio 
in the whole range of concentrations and interpolated by 
Eq. [2j whose coefficients a, /3, 7, 6 are given in Table |l] 
Coefficient 9, also given in Table |lj is the critical temper- 
ature of the miscibility gap. 

Homo-atomic pair interaction energies in Eq. [T] are 
given by the experimental cohesive energy of the pure 
elements, given in Tab. |l] according to Ecoh (j) = 
— J2i ^^^'^^^jj ■ The expressions of the hetero-atomic inter- 
actions (x, T) are then easily deduced from Eq.jljand 
[2j They are consequently simple parametric functions of 
the temperature T and local concentration x. This last 
quantity, x, around a pair including an atom on site i 
and an atom on another site j is naturally defined as 

(n) 

where p^f. — 1 when the kih site of the nth coordina- 
tion shell of site i is a Cr atom, and if it is a Fe atom 
or an empty site, i.e. a site outside the surface. The 
interaction range is restricted to second nearest neigh- 
bors with e-^^ = e^jV^, which has been found optimal. 
It is worth emphasizing that the resulting bulk phase di- 
agram is in very go od agreement with the most recent 
experimental review^^*^ While this model does not cap- 
ture the extraordinarily complex electronic structure of 
the bcc Fe-Cr alloys, it captures both the local nature 
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0.070 


-2.288 


4.439 


-2.480 


1400 
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4.10 



TABLE I. Parameters of the local-concentration and temper- 
ature dependent pair potential from Ref.'^. Experimental 
cohesive energies from Ref.'^ are given in eV per atom. 



of the interactions and the associated energetics, without 
empirical parameters. 

This letter focuses on the most stable surface of bcc 
iron, which has the orientation (lOOj^. It is modeled 
by a stack of 100 layers of 400 atoms each, in periodic 
boundary conditions. Interactions between periodic im- 
ages in direction (100) are prevented by a slab of vacuum. 
Special attention has been given to the choice of the size 
of the system in order (i) not to artificially hide the sur- 
face effects by a too large volume/surface ratio, and (ii) 
not to restrain the formation of precipitates by too small 
systems, ie to give the system the ability to precipitate 
and have the precipitates to interact with the surfaces. 
It induces that the bulk solubility limit near the surfaces 
can be slightly different from that of a pure bulk system. 
A perspective view of the supercell is shown in Fig. [l]a. 

As stated above, the model does not capture the effect 
of the surface on the electronic structure of atoms in its 
vicinity, such as expected stronger bonds. It captures the 
effect of the reduced coordination in terms of energetics: 
Eq.[3]implies that the more reduced the coordination, the 
more the local energetics are dependent on the remain- 
ing bonds. The various surface orientations only differ in 
the number of surface-induced dangling bonds, and thus 
in the strength of the surface effects described below. 
Conclusions are thus transferable to other orientations. 
Importantly, the effectiveness of our Hamiltonian allows 
to deal with a number of atoms that make it possible to 
finely tune the bulk and layer concentrations. Here, one 
atom accounts for the bulk concentration by less than 
10"^ at. % and the layer concentration by 5 • 10""^ at. %. 
This point is crucial as recent ab initio calculations have 
been limited to few layers and few atoms per layer, impos- 
ing large bulk and even larger planar concentrations. Our 
Hamiltonian is sampled by Monte Carlo simulations us- 
ing the Metropolis algorithm in the canonical and pseudo 
grand-canonical ensemble^^H xhe equilibrium state is 
considered reached after 10^ accepted permutations per 
site. 

We define the planar concentration Cp — 'J2f=i Ii l^p 
as the chromium content of each layer p parallel to the 
surface, with Np the number of atoms per layer (400 
here) and = 1 if site i d p contains a Cr atom, Qi — 
otherwise, p ranges from for the top surface layer to 99 
for the bottom surface layer. 

In figure [2j we plot the concentration profile Cp{p) of 
FeCr at 300 K (« 1/3 Tc) and various bulk concentra- 
tions Cb ranging from 0.02 to 0.98. Special attention is 
given to the temperature range that is of industrial and 
technological importance, i.e. for Cb below 0.3. The con- 



FIG. 1. Snapshots of the simulation cell containing the (100) 
FeCr surface at 300 K for bulk concentrations ct — 0.5 (left) 
and Ci, = 0.15 (right). The top surface layer is at the top of 
the figure. Fe atoms are shown in red and Cr atoms in blue. 
For Ci, = 0.15, only Cr atoms are shown. In both cases, the 
bulk phase separation a — a' occurs. 




FIG. 2. a) and b) Concentration profiles of a FeCr surface at 
300 K. Index p = indicates the top surface layer. Various 
bulk concentrations cj, are indicated: 0.02 (black); 0.05 (red); 
0.10 (green); 0.15 (blue); 0.30 (yellow); 0.50 (brown); 0.90 
(gray); 0.98 (violet); b) Only the first six layers are shown. 
Note the change in scale. 



centration profiles are highly non-linear functions of the 
bulk concentration. 

In order to get insight to these concentration profiles, 
semi-grand canonical Monte Carlo simulations have been 
performed, where the total number of sites and the dif- 
ference in chemical potential A/i between pure bcc-iron 
and pure bcc-chromium are kept fixed, while the bulk 
concentration is frec'^^. The evolution of the surface and 
bulk concentrations with respect to Afi at 300 K are plot- 
ted in Fig. [3] Three hysteresis loops are found, which are 
indicated in the figure by asterisks. They bring out three 
phase transitions. The first and stronger one, indicated 
by the black asterisk in Fig. [3]a., is an evidence of the 
well-known bulk phase-separation a — a' happening in 



FIG. 3. a) Evolution of the bulk concentration with the chem- 
ical potential A/i at 300 K. b) Evolution of the sum of the 
concentrations of the top surface and subsurface layers, ver- 
sus the same chemical potential. In both figures, the black 
asterisk indicates the bulk phase separation, the red and blue 
asterisks show the transitions in the subsurface and surface 
layers, respectively. 



bcc Fe-Cr alloys and discussed in the introduction. Note 
that the bulk solubility limit at low chromium concentra- 
tion is slightly affected by the presence of the surface. It 
causes the large variations in the density profiles in Fig. [2] 
more than 10 layers away from the surface for ct ~ 0.12. 
Snapshots of systems that undergo phase-separation at 
these concentrations are shown in Fig. [T] At higher dif- 
ference in chemical potential, as indicated in Fig. |3]b., 
two less visible phase transitions occur. Each transition 
is localized in a single layer. The first one corresponds to 
the subsurface layer transiting from pure Fe to pure Cr 
(indicated by the red asterisk in Fig. [sja., followed by an 
accompanying transition in the surface layer (blue aster- 
isk in Fig. |3^. and b.). The change in concentration of 
the two first layers is abrupt and discontinuous. It also 
gives insight to the emptiness in Cr of these layers. First, 
the difference in surface chemical potential of the two el- 
ements, which is proportional to the difference in surface 
energies of Fe and Cr, implies that Fe recover the layers 
where bonds are dangling. Indeed, surface energies of Fe 
are always lower than that of Cr for a given orientation. 
They range from 2.2 to 3.4 J/m^ and 3.2 to 4.2 J/m^ for 
iron and chromium, respectivelji^llEI!_ Secondly, and in 
relation, the chemical potential of surface atoms is much 
modified by the surface, which explains why the alloy do 
not phase separate at the same concentration than in the 
bulk. One could see them as two new alloying elements 
only present in the surfaces. 

As shown in Fig. [2jb and|4]a., the concentration of the 
third layer, C2, increases quickly with Cb, contrary to that 
of the first two layers discussed above: C2 gets from to 
0.2 when Cf, goes from to 0.1, which represents a relative 
increase in C2 of about 100 % at Ch ~ 0.1. Importantly, C2 
is at this point greater than the bulk solubility limit, so 
that phase separation would occur in the absence of the 
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FIG. 4. a) Relative evolution of the concentration of the third 
layer, C2, versus the bulk concentration, ct, at 300 K. b) Av- 
erage surface concentration, Cs — (cq + ci + C2 + C3) /4, as a 
function of ct, at 300 K (black squares), 450 K (red triangles) 
and 600 K (green stars). A black shaded line indicates the 
optimal concentration c^^^ ~ 0.12 at which ferritic FeCr steel 
is the most corrosion resistant. 



surface: its presence changes here the very nature of the 
Fe and Cr atoms, particularly their complex magnetic in- 
teractions as shown experimentallj^^l and theoreticalljff^, 
and consequently the perturbated alloy's thermodynam- 
ics. For higher bulk concentrations, C2 decreases sharply 
and becomes even depleted in Cr with respect to the bulk 
at Cb ^ 0.13. 

Finally, and importantly, we plot in Fig. |4]b the sur- 
face concentration, Cg, defined as the average concentra- 
tion of all layers in direct contact with the surface, i.e. 
Cs = Z]i=r Ci/«max with i^ax = 4 for Orientation (100), as 
a function of Cb at 300 K, 450 K and 600 K. Two regimes 
are clearly identified: (i) For bulk concentrations under 
0.12, the surface concentration increases with Cb, up to 
a narrow maximum between 0.09 and 0.12. As soon as 
Cb — 0.07, surface composition exceeds the bulk solubil- 
ity limit, where phase-separation would occur in the bulk 
alloy, (ii) For larger bulk compositions, there is a discon- 
tinuity in Cs, which is reduced to a flat regime Cg ~ 0.05, 
equivalent to that of a a — a' phase-separated bulk. In- 
deed, a' precipitation occurs as expected and discussed 
above to explain the large variations in the density pro- 



files of Fig. [2j and illustrated in Fig. [T] 
Between 300 K and 600 K, the profiles shown in Fig. |4] 
are almost temperature-independent. It reflects a sub- 
tle compensation between the temperature dependence 
of the order energy, more specifically of the energy asso- 
ciated with magnetism, and the entropic effects. It was 
identified in the bulk as the cause of the anomalously 
steep solubility limit of the Fe-Cr alloy at low temper- 

atureSj_and identified by Williams as an effect of niag- 
netisnPaiaa 

An atomistic explanation of the thermodynamic origin 
of the narrow optimum in corrosion resistance of stain- 
less steels emerges from the above results. The difference 
in surface chemical potential between Fe and Cr induces 
a strong Cr depletion in the first layers, where atoms 
have dangling bonds. This result can be understood as 
a surface effect resulting from the surface energies of Fe 
being always lower than that of Cr. The sub-surface lay- 
ers balance this local depletion by a strong enrichment 
in Cr, leading to three distinct regimes: (i) at low bulk 
concentrations, the ordering energy drives chromium into 
solution, far away from the surface: the resistance to cor- 
rosion is low and increases with Cr concentration, (ii) in 
a narrow range of bulk concentrations between 0.07 and 
0.12, enough Cr is present to strongly enrich the surface 
in average but not to exceed the a-a' solubility limit'^ 
the Cr content in the surface is maximum, (iii) at higher 
bulk concentrations, the bulk solubility limit is exceeded 
and most of available chromium bulk-precipitates in the 
a' phase, depleting the surfaces: stainlessness is lost. 

In light of the above results, in order to improve the 
stainlessness of the Fe-Cr system and its derivative al- 
loys, one may consider adding alloying elements that in- 
crease the solubility of Cr in Fe, without altering the 
anomalous thermodynamics of the Fe-Cr surfaces. The 
sharp decrease in the surface concentration would thus 
happen, nevertheless, but at higher bulk concentrations. 
Larger surface concentrations would be reached, induc- 
ing better protection against corrosion. One could also 
increase the quantity of chromium atoms at the surfaces 
without increasing the content in the bulk. This could for 
instance be done by the localized addition of chromium- 
rich nanoscale precipitates or dispersoids. Such strategy 
is already under investigation as it is also a promising 
route to strengthen these material d^ l ^^ l 
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